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Abstract 

The classical models used for describing the mechanical behavior of woven fabrics do not fully account 
for the whole set of phenomena that occur during the testing of such materials. This lack of precision is 
mainly due to the absence of energy terms related to the microstructural properties of the fabric and, in 
particular, to the bending stiffness of the yarns. The importance of the bending stiffness on the overall 
mechanical behavior of woven reinforcements, if already essential for the complete description of balanced 
fabrics, becomes even more important in the case of unbalanced ones. In this paper it is shown that the 
unbalance in the bending stiffnesses of the warp and weft yarns produces macroscopic effects that are 
extremely visible: we mention, for example, the asymmetric S-shape of a woven interlock subjected to a 
Bias Extension Test (BET). 

We propose to introduce a constrained micromorphic model and, simultaneously, a discrete model 
that are both able to account for i) the angle variation between warp and weft tows, ii) the unbalance in 
the bending stiffness of the yarns and iii) the relative slipping of the tows. 

The introduced constrained micromorphic model is rigorously framed in the spirit of the Principle of 
Virtual Work for the study of the equilibrium of continuum bodies. A suitable constraint is introduced in 
such micromorphic model by means of Lagrange multipliers in the strain energy density and the resulting 
constrained model is seen to tend to a particular second gradient one. The main advantage of using such 
constrained micromorphic model is that the kinematical and traction boundary conditions that can be 
imposed on some sub-portions of the boundary of the considered body take a natural and unique meaning. 

The discrete model is set up by opportunely interconnecting Euler-Bernoulli beams with different 
bending stiffnesses in the two directions by means of rotational and translational elastic springs. The 
main advantage of such discrete model is that the slipping of the tows is described in a rather realistic 
way. Suitable numerical simulations are presented for both the continuum and the discrete models and 
a comparison between the simulations and the experimental results is made showing a definitely good 
agreement. 


Introduction 

For decades textile composites made of woven fabrics have been successfully employed in aircraft and automo¬ 
bile engineering and they are gaining an even increased interest due to their excellent mechanical properties 
such as a very high specific-strength and excellent formability properties. Fibrous composite reinforcements 
present improved characteristics, namely high specific stiffness and strength, good deformability, dimensional 
stability, low thermal expansion, good corrosion resistance and many others. Among the quoted character¬ 
istics, the good deformability is what makes these materials perfect to be formed in various shapes with 
limited expenses. On the other hand, some complex behaviors of the woven fabrics, for instance the onset 
of wrinkling and slippage, limit the admissible deformations during the stamping operations and can render 
the modeling of such materials difficult to be achieved. Due to the complexity of the micro-macro behavior 
of fibrous composite reinforcements, the need of a comprehensive model for the prediction of the mechanical 
response of such materials during the forming represents a real scientific challenge (see e.g. [S]). 

The structure of the fabric is characterized by two main directions of woven tows (warp and weft) and, 
therefore, in those directions the fabric has a very high extensional rigidity. The way in which these fabrics 
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are weaved varies according to the different production methods and it is therefore possible to observe various 
schemes of weaving patterns such as those shown in figure for unbalanced fabrics. Each of these schemes 
leads to a different type of composite reinforcements and thus to different mechanical properties. Furthermore, 
the warp and weft of a fabric can be either balanced (with the same properties) or unbalanced (the warp 
and weft present different characteristics due, for example, to a different number of constituting fibers). Such 
unbalanced fabrics may be of use, for example, in all those engineering applications that require a material 
that has to be stressed in a main direction. In such situations the use of an unbalanced fabric conveys a real 
advantage in terms of material response and therefore an economic advantage as well. 



Figure 1: Schemes of weaving for unbalanced fibrous composite reinforcements 

In order to produce the fabrics, the yarns (that are themselves composed of thousands of small carbon 
fibers) are weaved together forming a very complex texture. In such a structure it is no wonder that the 
interaction between the yarns and their behavior at a mesoscale play a fundamental role on the overall 
response of the material. One of the main features that determines the properties of the fabric is the friction 
between the yarns that both prevents the slipping and generates the shear rigidity of the fabric. This shear 
stiffness fundamentally determines the behavior of the fabric being usually orders of magnitudes lower than 
the elongation stiffness of the yarns. The shear angle variation between warp and weft can hence be thought 
to be the principal deformation mode of fibrous composite reinforcements. However, the shear stiffness of the 
fabric and the elongation stiffness of the yarns are not enough to fully describe the response of such materials. 
The yarns, possess an in-plane bending stiffness that determines some peculiar behaviors of the material at 
the macroscopic scale (see e.g. HZllSl]). Particularly, in the case of an unbalanced fabric the difference in 
the thickness of the two families of yarns leads to peculiar responses. 

Notwithstanding the importance of a complete and effective modeling of fibrous composite reinforcements, 
the most common models for the description of the mechanical behavior of such materials fail to describe 
comprehensively their response. 

In this paper we propose to use a micromorphic model to describe the mechanical behavior of unbalanced 
fibrous composite reinforcements with an application to the Bias Extension Test, which is a very well known 
mechanical test in the field of composite materials manufacturing (see e.g. |lfll l4Tl l4^ I7n| b 

We try to limit at most the complexity of such model by introducing a unique kinematical scalar field ip 
in addition to the classical macroscopic displacement u. As usual in micromorphic models, the strain energy 
density of the considered material is supposed to depend on the first gradients of both u and (p. 

In order to describe at best the physics of the problem it is useful to interpret the variable p as the angle 
variation between warp and weft, so that a dependence of the strain energy density on Vtp actually accounts 
for the bending stiffness of the yarns. Indeed, it is possible to understand that when varying the angle from 
one value to a different one, this cannot be done too sharply, but a smooth gradient of the angle variation 
must occur that is physically associated to the bending of the tows. 

As a matter of fact, the angle variation between the two material directions mi and m 2 can be directly 
associated to the gradient of the displacement field (such angle variation is known to be the invariant ig = 
(mi, C ■ m 2 ), of the Cauchy-Green strain tensor C). In order to let the micromorphic variable p tend to the 
angle variation fg, we decide to use suitable Lagrange multipliers in the micromorphic strain energy density. 
The resulting constrained micromorphic model can be actually interpreted as a particular second gradient 
model since the first gradient of p tends to the gradient of ig which is itself a function of the gradient of u. 

As we will show in detail in the body of the paper, the proposed constrained micromorphic model is able 
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to account for 


• the basic mechanism of angle variation between warp and weft 

• different bending rigidities in the warp and weft directions 

• the decrimping and eventual relative slipping of the yarns. 

When applying such constrained micromorphic model to the description of the Bias Extension Test (BET) 
on an unbalanced fabric, we are able to recover the characteristic macroscopic asymmetric S-shape of the 
specimen (essentially due to the two different bending stiffnesses), as well as the deformation patterns of 
warp and weft tows at the mesoscopic scale. More precisely, for what concerns the slipping of the yarns 
which is experimentally observed, it is accounted for in our continuum model by means of the introduction 
of “equivalent elongations” in the two material directions. 

The reader may think that the fact of passing through a constrained micromorphic model to formulate 
a specific second gradient one is an artificial procedure that provides additional complexity to the modeling 
of the considered materials. Nevertheless, we will show that the fact of starting from a micromorphic model 
and to subsequently constrain it in order to tend to a second gradient one, is a more natural way to clarify 
the physical meaning of the considered problem. In particular, the fact of considering a micromorphic 
model avoids any confusion concerning the boundary conditions that can be imposed in the considered 
physical problem. We will see that we will be naturally led to impose kinematical boundary conditions that 
immediately take the precise physical meaning of imposing the displacement and the angle between warp and 
weft on given subsets of the boundary. Similarly, there will be a unique way to impose forces and double forces 
at the boundary of the considered micromorphic medium where tractions are assigned, and the introduced 
“force” will be seen to be directly related to the force which is measured by the employed testing machine. 

On the other hand, if we would have started directly from a second gradient model the way of imposing 
physical boundary conditions would have been much more complicated. In fact, it is known (see e.g. |52l 1561 
1^) that there is no a unique way to define forces and double forces in second gradient theories, but multiple 
combinations of such contact actions may be prescribed on the boundary which are all equally legitimate but 
which give rise to different boundary value problems. 

It is our hope to convince the reader that the physics of the boundary conditions which have to be imposed 
to model the BET is naturally suggested by the use of the used constrained micromorphic model. 

In a second time, we introduce a discrete model for the description of the Bias Extension Test on un¬ 
balanced fabrics. To do so, we use long Euler-Bernoulli beams with different bending stiffnesses in the two 
directions, which are suitably interconnected by rotational and translational springs in such a way that the 
main characteristics of the experimental Bias Extension Test are described. In particular, we are able to 
recover in a realistic way both the macroscopic asymmetric shape of the specimen as well as the pattern of 
the yarns at the mesoscopic scale which include shear strain, bending and slipping. Since the yarns constitut¬ 
ing the specimen never experience compression during the BET, we can sensibly affirm that the introduced 
model of interconnected Euler-Bernouilli beams actually describe, at least qualitatively, the overall behavior 
of the reinforcement. 

Even though the proposed discrete model only involve elastic elements (beams and springs), it is able to 
catch in a quite realistic way the main features of the deformation of unbalanced fabrics. A point that could 
be improved in this sense is to introduce some dispersion in the model (friction), following e.g. the methods 
presented in |361159| . in order to account for the phenomenon which is experimentally observed and which 
suggests that, after unloading, the specimen does not perfectly return to its initial undeformed shape. 

The results obtained via the discrete model allow us to comfort those obtained via the constrained 
micromorphic continuum model and they additionally allow to bring more light on the mechanisms of slipping 
that occur during the BET on unbalanced fabrics. 

The paper is organized as follows 

• In section 1 we set-up the equilibrium problem for both Cauchy and micromorphic continua. To do so, 
we introduce suitable spaces of configurations and of admissible variations whose structure is directly 
related to the kinematical boundary conditions that are imposed in the considered problem. We hence 
formulate the equilibrium problem for the considered (Cauchy and Micromorphic) continua by means of 
the use of the Principle of Virtual Work. Finally, we treat in more detail the problem of the equilibrium 
of a first gradient continuum by obtaining the irreducible form of the work of internal actions for such a 
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continuum. This irreducible form allows to individuate the type of boundary contact actions which can 
be introduced in Cauchy media, namely “forces” per unit area. Finally, we explicitly show a method to 
calculate external contact actions by means of the use of the Principle of Virtual Work: such method 
is based on a wise choice of particular test functions that allow to isolate forces and double forces on 
some specific parts of the boundary. 

• In section 2 we study the equilibrium of an unbalanced 2D fabric subjected to a Bias Extension Test. 
This study is carried out following three subsequent steps. First of all the physical problem (BET on 
unbalanced fabric) is precisely described with particular attention to the description of the imposed 
kinematical and boundary conditions (see e.g. [lOI [T^ 1551 Secondly, a discussion concerning the 
most appropriate continuum model which is needed to describe the physical phenomenon of interest 
is carried out. We propose a constrained micromorphic model to accomplish this task (see also [32]). 
Such model is able to account for i) the angle variation between warp and weft yarns, ii) the unbalanced 
bending stiffness of the two families of fibers and iii) the relative slipping of the yarns by means of the 
introduction of “equivalent elongations”. Last but not least the introduced constrained micromorphic 
model is seen to be able to introduce in a natural way the boundary conditions which are peculiar 
of the BET. Finally, we present the numerical simulations of the proposed constrained micromorphic 
continuum model to show that it is able to satisfactorily describe the experimental evidences. 

• In section 3 we introduce the discrete model for the description of the BET by means of the use 
of Euler-Bernoulli beams suitably interconnected by rotational and translational elastic springs. We 
show that the obtained results fit well the available experimental evidences both for what concerns the 
macroscopic and microscopic deformation patterns. In particular, the phenomenon of relative slipping 
of the yarns is unveiled in a rather precise way so allowing a wise interpretation of the equivalent 
elongations introduced in the continuum micromorphic model. 

We explicitly remark that section 1 contains an introductory theoretical treatise which is convenient to 
frame the considered mechanical problem in the framework of the Principle of Virtual Work. We believe that 
such a discussion is indeed very useful for those readers who want to make a clear connection between a neat 
mathematical formulation of the Principle of Virtual Work and the description of the mechanical phenomena 
that such principle can provide. Actually, an intelligible and self-consistent disquisition ranging from the 
mathematical setting to the mechanical interpretation of the Principle of Virtual Work is difficult to be 
found in academic articles, so that we believe that section 1 represents an added value for the present paper. 
Some of the introduced notations are voluntarily lightened with respect to those presented in classical more 
mathematical books (see e.g. |^) in order to make easier the connection between the mathematical and the 
mechanical aspects of the same problem. The reader who is not interested in establishing such connection 
can skip this opening section, directly passing to the mechanical setting-up of the considered problem. On 
the other hand, the reader who is interested in a more precise and general mathematical formulation of the 
problem can refer e.g. to |55| . 


1 The Principle of Virtual Work and the equilibrium of Cauchy and 
micromorphic continua 

In what follows we will call Cauchy continuum body a set of material particles occupying the volume B in 
its reference configuration and whose motion is described by means of a suitably regular, kinematical field 
u : B —>■ which we call the displacement field of the considered body. We denote bjQ 5u : B ^ TMf the 
virtual displacement field associated to the considered body. Finally, we will denote by dB the boundary of 
B. 

Generalizing this definition of Cauchy continuum, we will introduce a simple micromorphic continuum 
body by complementing the previously defined Cauchy continuum with a supplementary, suitably regular, 
scalar kinematical field tp : B —)■ M, that we generally call micro-motion. We denote bjj^ Sip ■. B ^ TM. the 
virtual variations of the kinematical field ip. We remark that, in the spirit of Mindlin |5bj and Eringen m, 
such supplementary kinematical field represents the motion of a microstructure which is embedded in the 

^Here and in the sequel we denote by TR^ the tangent bundle to the manifold R^. We recall that the virtual variation 6u of 
a displacement field u has the structure Su{X) = {u{X), (5ri7(JV)), where Sw{X) is a vector attached to the Eulerian point u{X). 

^Here and in the sequel we denote by TR the tangent bundle to the manifold R. We recall that the virtual variation S(p of a 
micro-strain field ip has the structure Sp{X) = {p{X),S'il)(X)), where 5'0(V) is a vector attached to the image point p{X). 
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considered body. Such micro-motion is, in principle, completely independent of the macroscopic motion of 
the matrix. Nevertheless, in some cases of physical interest, as the one which will be analyzed in this paper, 
it is worth to relate such micro-descriptors to the derivatives of the macroscopic displacement field. As we 
will show later on in much more detail, this can be done by constraining the introduced micromorphic model 
with suitable Lagrange multipliers to be added in the strain energy density, or equivalently, by introducing 
a penalty term in the energy itself. We want to stress the fact that we prefer to keep a more general 
micromorphic model and to subsequently constrain its strain energy density in order to let it tend to a 
second gradient one. This choice is preferable if one wants to interpret in a unique way the external actions 
of the considered continuum. In fact, (see e.g. |52l IMl I57| l if one starts directly from a second gradient 
energy, the interpretation of boundary contact actions, namely forces and double-forces, depends on the type 
of manipulation which is done on the work of internal actions by means of procedures of integration by parts. 
More particularly, if one decides to stop at a given level of integration by parts, or to continue further to the 
subsequent level, the definition of force and double force is not the sam^ When considering micromorphic 
continua in which only first gradient of the introduced kinematical fields appear in the strain energy density, 
only one level of integration by parts is possible and then the boundary contact actions are uniquely defined 
and take immediate physical meaning when framed in the considered physical problem. 

The equilibrium of a Cauchy continuum body subjected to given boundary conditions can be studied 
by means of the Principle of Virtual Work. Such fundamental principle of Mechanics states that a body, 
subjected to specific external actions, is in equilibrium if the work of internal actions is balanced by the work 
of external actions. In formulas, we say that a displacement field u* is an equilibrium configuration if 

{u *, Su) + {u *, Su) =0, (1) 


for any compatible 6u. In most cases, the external and internal works can be seen as the first variation 
of suitable functionals A™* (u) : Q M and (u) : Q —)■ K, so that the Principle of Virtual Work Q 
actually implies the minimization of a functional A := -I- More specifically, we can write 

{u, 6u) + (u, 5u) = SA{u, Su) := lim ^ ^ ( 2 ) 

(->■0+ t 

where we denoted by 8A the first variation of the functional A, where u G Q, Su € and the sets Q and T„ 
will be defined in more detail later on. 


Suitably generalizing the definitions given above for first gradient continua, the Principle of Virtual 
Work can be reformulated for the introduced micromorphic continuum by saying that a couple (u*, ip*) is of 
equilibrium if 

jyint (^y*^ ip*^ (u*,(f*,Su, Sip) = 0, (3) 

for any compatible {Su,Sp>)- 

Again, the external and internal works can be seen as the first variation of suitable functionals A*"* {u, ip) : 
Q X D ^ M. and {u,ip) : Q x D ^ M., so that the Principle of Virtual Work ([^ actually implies the 
minimization of a functional A := -I- More specifically, we can write 


P™* {u,ipAu, Sip)+V^^* 


(u, ip, Su, Sip) = SA(u, ip, Su, Sip) 


A(u + tSu,ip + tSip) — A(u,ip) 

lim - 

(->■0+ t 


(4) 


where we denoted again by SA the first variation of the functional A, where (u, ip) € Q x D, (Su, Sip) S x 
and the sets Q, D, r„ and T^p will be defined in more detail later on. 

Suitable generalizations of the Principle of Virtual Work can be introduced in order to account for inertial 
and dissipative effects, but, since we deal in this paper only with static problems, we refrain here to present 
such more complex framework. 

The most fundamental questions which have to be confronted to properly set up a mechanical theory by 
means of the Principle of Virtual Work is to establish: 

3 

We need to mention the fact that no common agreement is currently available concerning the choice of different but equally 
legitimate sets of boundary conditions deriving to different levels of integration by parts. Nevertheless, this is what is found e.g. 
in |52l 1561 [57| and we tend to adopt this viewpoint in the recent times. 


“^We remark that, depending on the conventions which are used for the signs in the definition of the work of internal and 
external actions, slightly different versions of the Principle of Virtual Work can be found in the literature. 
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• the constitutive form of the work of internal actions in terms of the displacement and, eventually, of 
the micro-descriptor (such constitutive choice is related to the intrinsic nature of the medium that one 
wants to study), 

• the expression of the work of external actions, which allows to establish how the external world acts 
on the considered medium and to define the concept of force, double force, or other more complex 
interactions. 

As a matter of fact, we can imagine to act on the boundary of the considered body by imposing either 

• kinematical (or essential or geometric) boundary conditions: the displacement and/or eventually the 
micro-descriptor are assigned on some portion Hk of the boundary dB, 

• traction (or natural) boundary conditions: forces and/or, eventually, other more complex external 
interactions are assigned on some portion S-p of the boundary dB. 

In order to be more general, we can introduce 

• The surface ^Ki on which displacement is assigned and the surface on which we can eventually 

assign the micro-descriptor ip. Clearly, IJ E/^^j the two sets Ejf^ and Ejf^ may eventually 

partially or totally overlap depending on the considered physical problem. 

• The surfaces E^i on which we assign forces and E^^ on which we can assign more complex interactions 
and that can as well partially or totally overlap depending on the physical problem in study. We also 
have Et = Et, U ^T2- 

If, in the considered physical problems (this may eventually arrive in micromorphic theories), mixed 
kinematical-traction boundary conditions are applied, also the sets and Ey^ and/or Ejf^ and E^i may 
eventually have non-vanishing intersection. 

We also explicitly mention that, in order to be able to recover all the possible external interactions, each 
of the introduced sets ^Xi, '^k^i ^Ti and Ey^ niay collapse in the empty set or can cover the whole boundary 
dB. Finally, we remark that, by definition, we set Hx (J E^ = dB. Finally, we also mention that we will 
neglect external volume forces in the following treatment. 

1.1 Space of configurations and spaces of admissible variations 

Depending on the intrinsic nature of the considered body (first gradient or micromorphic), the expressions 
for the work of internal and external actions take specific forms which will be better specified later on. 

Independently of the specific form taken by the internal and external work, we want to underline here 
that the problem of finding the equilibrium configuration of a given continuum (Cauchy or micromorphic) 
subjected to specific boundary conditions reduces to the problem of finding, in suitable sets, the kinematical 
fields which satisfy the Principle of Virtual Work (([^ or (§) for any virtual admissible variations. 

We start by introducing the equilibrium problem for a Cauchy continuum, defining a suitable set Q, 
called space of configurations of the considered medium which contains information about the kinematical 
constraints which must be verified by the displacement field. More precisely, we define the space of configu¬ 
rations for a Cauchy continuum a^ 


Q = {u I u = u on Eifj}, (5) 

where u is a suitably assigned function. Roughly speaking, the set Q represents the set in which we look 
for the solution of our minimization problem and contains only those displacement fields which satisfy the 
imposed kinematical boundary conditions. 

On the other hand, we define the set of admissible variations as 

T„ = {Su \u-\- 5u ^ Q} . (6) 

®We remark that the set Q should also contain informations concerning the desired regularity on u, but we limit ourselves 
here to talk about “suitably regular” functions. 
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We explicitly remark that, being u = u on in order to have 6u belonging to the set of admissible 

variations, one must have that u + Su = u on ■ This clearly implies that = 0 on TiKi > and hence the 
set of admissible variations takes the form Tu = | = 0 on S/fjj.With the introduced notations, we 

can formulate the equilibrium problem for a Cauchy continuum as: 

Find u* G Q such that {u*, Su) + 7^®^*(it*, Su) = 0 for any Su € T„.. 

The setting-up of the equilibrium problem for a micromorphic continuum can be obtained, suitably 
generalizing what done for Cauchy continua. In particular, a supplementary set £> must be defined to specify 
the space of configurations for such continuun:]^ 

D = {(fil ip = (f> on EicJ, 

where is a suitably assigned function. Roughly speaking, the set D represents the set in which we look for 
the solution for the micro-motion of our minimization problem and contains only those micro-motions which 
satisfy the imposed kinematical boundary conditions. 

Moreover, we define a supplementary set of admissible variations as 

= {Sip \ p + Sp & D} . 

With the introduced notations, we can formulate the equilibrium problem for a micromorphic continuum as: 
Find {u*,p*) G Q X D such that V™*{u*,p*,Su,Sp) + P^^*^{u*,p*,Su,Sp) = 0 for any {Su,Sp) G 

X . 


1.2 The example of the equilibrium of a first gradient continuum 

As it has been previously pointed out, in order to establish the equilibrium problem for a given continuum 
body subjected to specific external interactions, the expressions of both the work of internal and external 
actions must be specified. 

For a first gradient continuum the work of internal actions is defined through the definition of the action 
functional A™* which can be introduced in the static case a^ 

= - [ W (Vu) dv, 

J B 

where W is the strain energy density which, in a first gradient model, constitutively depend only on the first 
gradient of displacement. 

In the case of first gradient theories we can hence recognize that the work of internal actions can be 
written as|3 




SW(Vu) dv 




dv 







n 



ds, 


(7) 


where to obtain the last identity the divergence theorem has been used. Equation Q furnishes the irreducible 
expression of the work of internal actions for a first gradient continuum. This means that no more integrations 
by parts can be performed to ulteriorly manipulate this expression of 7^™*. It can be remarked from such 
irreducible form of the internal work that, as far as the boundary dB is concerned, only quantities expending 
work on the virtual displacement Su (i.e. forces) can be recognized. It is for this reason that, based on the 
validity of the Principle of Virtual Work, we can affirm that the only boundary external actions that can 
be sustained by a first gradient continuum are forces per unit area, i.e. external actions expending work on 

®We remark that the set D should also contain informations concerning the desired regularity on ip, but we limit ourselves 
here to talk about “suitably regular” functions. 

^Classically, in the dynamic case the internal action functional is defined as the space-time integral of the Lagrangian density 
C = T — W, where T is the kinetic energy density. In the particular case of statics, the minus sign remains after the due 
simplifications. 

®The operator Div stands for the classical divergence operator. Being A a tensor field of any order n > 0, we define its 
divergence as the n — 1 tensor (DivA)^^ ^ ^ = A^^. Finally (a, 6) = indicates the scalar product 

between two tensors of any order n > 1 and the Einstein convention of sum over repeated indices is used. 
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Su. These observations are at the origin of the introduction of the work of external actions for first gradient 
continua in the forir0 

V^^\u,5u)=l {f,Su)ds, (8) 

where / : i? —>■ is a suitable function assigned on . 

Once assigned the specific form for the work of internal and external actions, the equilibrium problem for 
a first gradient continuum can be hence reformulated as follows: 

Find u* G Q such that Ju) + 7^®^*(rt*, 5u) = 0 for any Su G , 

where now 7?™*, pext^ q given by Q, (|^, ([^ and (|^ respectively. 


1.2.1 Evaluation of the reaction force corresponding to an imposed boundary displacement. 


Once that the solution u* of the equilibrium problem for a first gradient continuum has been found following 
the steps presented in subsection |1.2[ it may be interesting to know which is the reaction force acting on 
Exi that balances the displacement u which has been imposed on the surface itself. In other words, 
we are looking for the force that one should apply on the portion of the boundary in order to produce 
the displacement u on and the displacement field u* within the considered body. In order to answer to 
this question, it is possible to pass again through the Principle of Virtual Work, but imagining now that a 
work of external forces must be introduced also on the portion of the boundary Yki ■ We hence re-define the 
external work ([^ by adding an additional term as follows 

= / {f,5u)ds+ f {fR,6u)ds. (9) 


The reaction force can hence be calculated by writing the Principle of Virtual Work as 


Div 


/ dW 


Su) dv — 


ldB=Y.T^ USjfi 


! dW 


Sii\ ds + f {f,Su)ds+ f {fii,Su)ds = 0 

/ J Sti J 


( 10 ) 

and imposing that it must be valid for any arbitrary displacement Su. Such identity must be satisfied, in 
particular, for a virtual displacement field Su which is such that 


• Su is constant on Yki 

• Su is continuous on dB 


• Su is vanishing outside except on a suitably small region in order to preserve continuity with the 
imposed displacement at the boundary. 


For such particular test function, the Principle of Virtual Work (101 implies 


The constant test field Su can hence be simplified on the two sides and, introducing the quantity R := 
/s Ir ds, we can finally write 

f dW 

We call R the reaction force associated to the imposed displacement u on ■ We will see that the fact of 
computing the reaction forces on the basis of the procedure shown here may be of interest for comparing the 
results of the performed numerical simulations to the experimental data. In fact, when a standard testing 
machine measures a “force” associated to an imposed displacement, it is actually measuring a resultant force 
on the considered boundary. 

We explicitly remark that, an equivalent way to calculate the reaction force can be found considering the 
work of internal forces before integrating by parts as given in the first equality in expression 0,i.e. 


^We recall once again that, to the sake of conciseness, we suppose in this paper that bulk external actions are vanishing. 









With reasoning analogous to the previous ones, if we write the Principle of Virtual Work evaluated in the 
equilibrium solution u*, we have 

-J + J {fR,6u)ds = 0, 

which must be satisfied for any Su. If we hence choose a particular test function Su such that 


1 . 6u is constant on 

2. 6u is continuous on dB 

3. 6u is an arbitrarily assigned, non-vanishing function outside ^Ki 


we can evaluate the reaction force by noticing that 


/ dW 


=/a aw 


Vi5m ) dv — 


if, 5u) ds. 


( 12 ) 


We explicitly remark that Eq. (12 1 is a scalar equation, which means that the three components of the vector 


R cannot be directly calculated only using such equation. Nevertheless, if we choose suitable test functions Su 
which are aligned with one of the directions of the used reference system, we can arrive to evaluate separately 
the components of R. In particular, let 23 }? be an orthonormal basis with respect to which we want 

to evaluate the three components Ri = (i?, e^) of the reaction force R. We choose a test function Su which 
possesses all the characteristics previously listed except that the point 1 . is replaced by 


Su is equal to Ci on 'Eki 


In this case the component Ri can be easily calculated according to Eq. (12). 

Clearly, if the problem is well formulated, the value of the reaction R calculated with expression (111 
must coincide with the one calculated using expression (12). The second way of evaluating the reaction 
force has many advantages, especially when we want to do it numerically. Indeed, it is always more stable 
to numerically evaluate a volume integral than a surface integral. It is for this reason that, for the case of 
micromorphic continua shown in the next section, we will limit ourselves to show how to calculate reaction 
forces and double-forces, by passing trough the evaluation of volume integrals. 


2 Equilibrium of a 2D unbalanced woven fabric modeled as an or¬ 
thotropic micromorphic continuum 

It has been known since the pioneering works by Piola [72], Cosserat [Tl], Midlin [^, Toupin |5^, Eringen 
m, Green and Rivlin m and Germain cm EEI that many microstructure-related effects in mechanical 
systems can be still modeled by means of continuum theories. It is known since then that, when needed, 
the placement function must be complemented by additional kinematic descriptors, called sometimes micro- 
structural fields. More recently, these generalized continuum theories have been widely developed (see e.g. 
iiHiiiniimiiiiiiiiiiisaisniiMiiMiiTniisn]) to describe the mechanical behavior of many complex systems, 
such as e.g. porous media Uni Eni IH2 ES], capillary fluids dn Hi na El EH, exotic media obtained by 
homogenization of heterogeneous media IS EH I84| . Interesting applications on wave propagation in such 
generalized media has also gained attention in the recent years for the possible application of this kind of 
materials to passive control of vibrations and stealth technology (see e.g. 

In this section, we will frame the problem of determining the equilibrium of a 2D, unbalanced, woven 
reinforcement modeled as a micromorphic material, in the spirit of section With this precise aim in 
mind, we will first present the particular physical problem that we want to study, then we will propose 
suitable expressions for the work of internal and external actions and we will introduce the adapted spaces 
of configurations and of admissible variations. We will finally formulate the equilibrium problem for the 
introduced micromorphic continuum and we will solve it by means of suitable numerical simulations. 
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2.1 The physical problem: bias extension test on an unbalanced fabric 

We consider in this subsection the description of the physical problem that we want to analyze and for which 
we will develop a suitable mechanical second constrained micromorphic model in the next subsection. The 
material that we consider here is an unbalanced fabric, i.e. a fibrous reinforcement in which the warp and 
weft yarns have very different thickness due to the fact that they are constituted by a significantly different 
number of fibers. This unbalance is reflected in a different mechanical behavior of the two orders of yarns. In 
order to sketch a simplified description of the mechanical behavior of unbalanced fabrics we can notice that 

• the yarns can be supposed to be quasi-inextensible in both the warp and weft directions, due to the 
very high tensile resistance of the constituent carbon fibers. This means that, notwithstanding the 
different thickness of the yarns, we can suppose that they do not sensibly elongate. Actually, an initial 
apparent elongation of the yarns due to decrimping can be eventually observed when testing the woven 
material, but it is reasonable to suppose that such apparent elongation, if present, has eventually the 
same characteristics in both directions. 

• The warp and weft being strongly unbalanced, we can infer that they possess very different bending 
stiffnesses. More particularly, the thick yarns are sensible to exhibit a higher resistance to bending 
than the thin ones. 

• Depending on the nature of the externally applied load and/or boundary conditions, the yarns can 
experience some relative slipping. More precisely, it is possible that the contact point between two 
yarns can move during the deformation of the macroscopic piece. 

We are interested in this paper to the description of a Bias Extension Test on an unbalanced fabric of 
the type described above. More precisely, the specimens object of the study are rectangular carbon fiber 
interlocks (the height must at least 2.5 times longer than the basis) with unbalanced yarns oriented at ±45° 
with respect to the long side of the specimen (see Fig. |^. The two shorter edges of the specimens are 
clamped in suitable devices which assure the following kinematical boundary conditions 

• vanishing displacement on Ei, 

• constant assigned displacement uq on E 2 , 

• angle between warp and weft blocked at 45° on Ei and E 2 (i.e. vanishing angle variation between the 
two families of yarns during the motion of the fabric). 



Figure 2: Experimental set-up for a bias extension test on an unbalanced fabric: it must be IE 4 I > 2.5 |Ei|. 
For the considered experimental test we have |Ei| = 70 mm and IE 4 I = 220 mm. Moreover, the specimen is 
15mm thick, but we suppose that this has no influence on the results, i.e. no displacement neither deformation 
occur out of the plane of the fabric. 

The experimental result of the performed mechanical test is a S-shaped macroscopic deformation as the 
one presented in figure 

The obtained asymmetric shape is related to the fact that the material properties and, in particular, the 
bending stiffnesses of the yarns are not the same in the two privileged material directions. Moreover, due 
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Figure 3: Experimental deformed shape for an imposed displacement uq = 56 mm. 


to the applied boundary conditions, some non-negligible slipping of the yarns is also observed, above all for 
what concerns the central part of the specimen in which yarns with two free ends are located. 

We will more precisely describe later which is the precise pattern of the yarns inside the considered 
specimen. On the other hand, we want to point out here that, indeed, a sensible differential bending can 
be observed in the considered specimen. In fact, the thin yarns are seen to be sensibly bent in some thin 
transition layers, while the thick yarns do not bend at all as highlighted in Fig. Such differential bending 
is certainly at the origin of the asymmetric deformed shape. 

2.2 The equilibrium problem for an orthotropic micromorphic continuum sub¬ 
jected to Bias Extension Test 

In this subsection we want to set up a consistent continuum mechanical framework which is able to account 
for 

• the presence of two preferred material directions inside the macroscopic body, 

• the unbalance of the fabric which, for what said in subsection [2d] basically means that the two families 
of yarns have different bending stiffnesses, 

• the relative slipping of the yarns. 

In order to do so, we need to remark that the desired continuum model (constrained micromorphic) must 
be such that 

• the chosen first gradient constitutive laws account for the orthotropy of the considered medium, 

• it allows to describe the curvature of material lines. Moreover, the resistance to curvature must be 
different for the two family of yarns, which means that the introduced constitutive laws must account for 
different bending stiffnesses for the two families of yarns. In order to reduce the adopted micromorphic 
model to a second gradient one, suitable constrains must be introduced in order to let higher order 
derivatives appear in the strain energy density, 

• notwithstanding the quasi-inextensibility of the carbon yarns, it includes equivalent elongation modes in 
the yarns directions in order to indirectly account for the slipping of the fibers. The fact of considering 
such elongation modes, also allows the model the possibility of describing the initial decrimping of the 
yarns. 


In order to formulate the equilibrium problem for the considered continuum subjected to a Bias Extension 
Test in the general framework introduced in section]^ we need to specify the particular form that 7?“*, 

Q , D, Tu and T^p take in the treated example. 
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Figure 4: Differential bending of the thick (b) and thin (c) yarns as observed in the experimental deformed 
shape (a). 


2.2.1 Work of internal actions for a 2D, orthotropic, constrained micromorphic continuum 

We start by specifying the expression of the work of internal actions which is suitable to describe the 
deformation of the considered system. To do so, we start by recalling classical results for first gradi¬ 
ent, hyperelastic orthotropic continua which prescribe the functional dependence that the strain energy 
density of an orthotropic continuum must have on the Cauchy-Green strain tensor C = ■ F (see e.g. 

la 0 [la ig ig Ell Ea [Ml Ell iig i87j) . 

Here and in the sequel, we denote hy x ■ B ^ the placement function associated to the considered 
body that associates to any material particle X G B its current position x in the deformed configuration and 
which can be related to the displacement field by means of the relation x = u + X. Moreover, we denote by 
F = Xx the space gradient of the introduced placement field. 

We start by assuming that the strain energy density can be given in the form 

W(C,X(p) = Wi(C) + Wii(Xip). (13) 

Representation theorems for 3D orthotropic first gradient materials are available in the literature (see e.g. 
Eiiigiini), which state that the functional dependence of the strain energy density on the strain tensor 
C must be given in terms of its invariants io ■= {*!)* 4 j*6)*8)* 9 )*io}) where the introduced invariants are 
defined in table ^ in which also their physical interpretation can be found. 

Explicit expressions for the strain energy potential as function of the invariants io which are suitable 
to describe the real behavior of orthotropic hyperelastic materials are difficult to be found in the literature. 
Certain constitutive models are for instance presented in |46| . where some polyconvex energies for orthotropic 
materials are proposed to describe the deformation of rubbers in uniaxial tests. Explicit anisotropic hypere¬ 
lastic potentials for soft biological tissues are also proposed in |M] and reconsidered in I3IHI] in which their 
polyconvex approximations are derived. Other examples of polyconvex energies for anisotropic solids are 
given in [SS] • It is even more difficult to find in the literature reliable constitutive models for the description 
of the real behavior of fibrous composite reinforcements at finite strains but some attempts can be for instance 


12 















recovered in Hill]. Furthermore, the mechanical behavior of composite preforms with rigid organic matrix 
(see e.g. |2811551 [Ml I69| l is quite different from the behavior of the sole fibrous reinforcements (see e.g. [T3]) 
rendering the mechanical characterization of such materials a major scientific and technological issue. 


Invariant 

Expression 

Meaning in terms of deformation 

ii 

tr(C) 

Averaged changes of length 

H 

mi ■ C ■ mi 

Local stretch in the direction mi 


mg - C ■ m2 

Local stretch in the direction mg 

is 

mi ■ C ■ mg 

Angle variation between the directions {mi, mg) 

*9 

mi ■ C • m3 

Angle variation between the directions {mi, m2) 

*10 

mg - C ■ m2, 

Angle variation between the directions {mg, m2) 


Table 1: Invariants of the Green-Lagrange strain tensor in the orthotropic case. The vectors mi and m 2 are 
unit vectors in the two privileged directions of the material and m 3 := mi x m 2 . 


For the particular 2D case that we study here we assume that the first gradient energy takes the following 
particular constitutive form 


Wi{C) = [iVh - 1)' + (v^ - 1)"] + ^Kshil 


(14) 


where the parameters K^i and Kgh are assumed to be constant. The constitutive choice (14) can be motivated 


by noticing that the invariants ig and iig represent out-of-plain angle variations of the yarns which are not 
considered in the present 2D case. Moreover, in a first instance, the isotropic invariant ii can be considered 
to be uninfiuential in the considered orthotropic case. 

It must be remarked that this simple quadratic choice for the first gradient strain energy density, even 
if providing geometric non-linearities, could be not sufficiently general to describe larger deformations for 
which more complex hyperelastic constitutive laws should be introduced. More than that, since for very 
large strains the integrity of the material starts to be affected due to the excessive slipping, there is no 
interest in attempting the modeling of the targeted unbalanced materials after a given strain threshold. In 
order to precisely identify the maximum strain that the material can withstand before failure due to excessive 
slipping, more experimental campaigns should be carried out. We limit ourselves here to remark that: 


• After a first threshold the material behavior starts to present a softening (see experimental pattern in 
Fig. § which can be directly related to slipping. To model such behavior, more general hyperelastic 
laws with respect to the one presented in this paper should be introduced, but this falls outside the 
scope of the present work. 

• If the experiment is prolonged, the slipping becomes so important that the integrity of the material 
starts to be affected and some yarns are pulled out of the specimen. In this case, both discrete and 
continuum model loose their predictability. 


In summary, we are saying that, with the considered constitutive choice and when remaining in the mod¬ 
erate strain regime, the main first gradient deformation modes allowed in the deformation of the considered 
material are 


• the angle variation ig between the warp and weft direction 

• the equivalent elongations 14 and ig in the directions of the warp and weft which account for decrimping 
and, eventually for slipping. 


As for the micromorphic energy, we make the following particular constitutive choice 


Wu{V^) = \{a,V^)^ 


- + a2(p,2 ), 


(15) 


where the vector a = ( 01 , 02 ) is the vector of constant micromorphic elastic parameters whose components 
oi and 02 have to be different to account for the unbalance of the microscopic characteristics of the material. 
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Moreover, we denoted by Lp^i and ip ^2 the space derivatives of Lp with respect to the space coordinates in the 
directions mi and m 2 (warp and weft) respectively. 


We explicitly remark at this point that we want to give a precise physical meaning to the micro-descriptor 
Lp in order to catch at best the experimental behaviors described in subsection |2.1| We have pointed out 
that a supplementary deformation mode with respect to classical first gradient ones must be introduced in 
order to catch all the physics of the considered problem. More particularly, additionally to angle variations 
(ig) and equivalent elongations in the yarns’ directions (14 and ig), we need to account for the bending of the 
two families of fibers. It is known (see also miisg) that bending strains of the fibers inside the considered 
macroscopic specimen can be accounted for by introducing second derivatives of the displacement field in 
the strain energy density. In particular, such bending of the yarns can be related to the space gradient of 
the angle variation ig: if sharp variations of angle between warp and weft occur in small regions inside the 
specimen, it means that the yarns must necessarily bend in order to rapidly change their direction and give 
rise to such angle variation. 

In the light of such remarks, it is sensible to suppose that the introduced micro-descriptor Lp must be 
indeed related to the angle variation fg: if, for example, we let ip tend to ig, then expression (15 \ for the strain 
energy density accounts for space derivatives of the angle variation and hence, finally, for the bending of the 
two families of yarns. Indeed, (see also |17[ I32p . it may be rather easily inferred how fg_i can be interpreted 
as the bending of the yarns disposed in the mi direction and, analogously, fg _2 represents the bending of the 
yarns aligned in the m 2 direction. 

Based on the physics of the problem discussed in the previous subsection, we do not introduce second 
gradient effects related to the gradients of the other invariants. We are then excluding that sharp spacial 
changes of elongation occur in the considered material. 


At this point, the reader may believe that the fact of considering a micromorphic medium is redundant 
for treating the considered problem of the Bias Extension Test, since a second gradient model could have 
been directly introduced, instead of constraining a micromorphic model to become a second gradient one. 
Nevertheless, the intermediary step of passing through a micromorphic model is essential, at least for two 
reasons 


• the imposed boundary conditions take a precise and unique meaning 

• the numerical implementation of the considered problem is more easily treatable since lower order 
differential equations are involved. 


The first point of the unique meaning of the imposed boundary conditions is crucial if one wants to deal 
with a model which has an easily recognizable physically grounded interpretation. In fact, as far as second 
gradient theories are concerned, the boundary conditions that can be imposed may take different, but equally 
legitimate, forms for the same physical problem (see e.g. [52]): for example, in a second gradient theory, a 
given angle can be imposed either by directly assigning the angle or by suitably choosing the components of 
the normal derivative of displacement on the boundary. Depending on whether one choice of the kinematical 
conditions or the other one is made, the dual traction counterparts (dual of the angle variation or of the 
normal derivative) have different expressions and the definition of the force can be also shown to be non¬ 
equivalent in the two cases. Such non-uniqueness of the way of imposing second gradient boundary conditions 
is directly related to the number of integration by parts which one decides to make in the expression of the 
internal work: in second gradient theories, the second gradient of the virtual displacement can be integrated 
by parts twice, by making use of the standard divergence theorem and of the surface divergence theorem. 

On the other hand, micromorphic models only involve first gradients of the introduced kinematical fields, 
so that only one level of integration by parts can be conceived (only the standard divergence theorem is used 
in a micromorphic model). This fact, avoids any sort of indeterminacy for the imposable boundary conditions 
when micromorphic models are considered (see e.g. also 0 ). 


In order to implement the fact that we want to constrain the micromorphic energy (131 to a second 


gradient one, we need to introduce a suitable Lagrange multiplier A with associated strain energy density 

VEa 

WA{A,is,ip) = A{ip - is) (16) 


In this constrained framework the global energy of the system is not simply (13), but must be comple¬ 


mented with this additional coupling term and in the considered 2D case becomes 


W{i4,i(i,is,ip,A^(p,A) = Wi{i4,ie,is) +Wii{W(p) -k Wa(A, jg, v?)- 


(17) 
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We explicitly mention that, the introduction of the Lagrange multiplier A does not modify the definition of 
the set D. 


According to the constitutive choices (14) and (16), the work of internal actions of the considered 2D, 


orthotropic, constrained micromorphic continuum can be written as 

= -^ [ {Wi + Wh + WK)dv = - [ {SWi + 6Wn + SWa) dv. 

Jb Jb 

Computing the first variation of the introduced internal action functional, we can write 

=~ Sis + A^^SVip + -A^Sifi + -a^SA 

Jb I ^^4 ^^6 \ ^^8 ^'^8 J U\<^ o(p uA 


dv. 


which for the particular constitutive choice made to treat the physical example considered in this paper, 
takes the particular form 


^int _ 

'II — 


^1 - Si4 - ^Kei ^1 - -A j Sie - {Kshh - A) Sis 


(18) 


— {a,W(p) {a,WS(p) — AS(p — {(fi — is)SA 


dv. 


Considering the definitions of the invariants given in table (Q, it can be checked that 
Si4 = mi ■ S (F^ • F) • mi = 2 (F • mi) • {SF ■ mi) 


Sis = m2 • 5 (F^ • F) • m2 = 2 (F • m 2 ) • {SF ■ m 2 ) 


Sis = mi ■ S (F^ • F) • m2 = (F • m 2 ) ■ {SF ■ mi) + (F • mi) • {SF ■ m 2 ) , 

so that, replacing such expressions in the first three terms of the internal work ( [I^ , recalling that SF = 
^('^x) = V(<5ri), suitably integrating by parts each term and using the divergence theorem we finally get 


= ( Div 


Kei ( 1 - ) (F • mi (g) mi) 


Div 


Kei ( 1 - ) (F • m 2 (g m 2 ) 


, Su) dv 


+ / ( Div [ {Ksh is — A) F ■ (mi g) m 2 + m 2 g) mi) ], Su) dv — {ip — is)SA dv 

J B J B 


IdB 


Kei ( 1 - ) {mi ■ n) + {Ksh is - A) {m 2 ■ n) 


(F • mi) , Su) ds 


(19) 


IdB 


Kei ( 1 - ) {m 2 ■ n) + {Ksh is - A) (mi • n) 


(F • m 2 ) , Su) ds 


+ / [a • V (( a, )) - A] Sip- {a,Vip){a,n) Sip 
J B JdB 


where we denoted by n the unit normal to the Lagrangian boundary dB. We explicitly remark that imposing 
arbitrary variations SA implies the desired constraint ip = is- 

The expression ( |19[ ) is the irreducible form of the work of internal actions for the considered constrained 
micromorphic medium. We explicitly remark that, the boundary contact actions intervening in the considered 
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problem expend work on 5u (forces) and on 5ip (double-forces). As far as the boundary dB is concerned, it 
can be inferred from Eq. that the dual quantity to the virtual displacement 6u contains informations 
about the first gradient deformations ^ 4 , iq, is as well as the Lagrange multiplier A. We are in some way 
saying that the internal force which would eventually balance an externally applied force involve deformation 
mechanisms related to elongations of the warp and weft, angle variations, but also (through the Lagrange 
multiplier A) some microstructure-related deformation modes. On the other hand, the internal double-force 
(couple) which would eventually balance an externally applied double force, involve deformation mechanisms 
associated to the local bending of the yarns (first derivatives of the angle variation ip = is). 

Analogously to what done in subsection |1.2.1[ we could use this irreducible expression of the work of 
internal forces to suitably calculate the reaction forces and double-forces on the part of the boundary where 
we assign the kinematical constraints. Nevertheless, since in our numerical simulations we only use the 
second method to calculate such reactions (passing through the evaluation of volume integrals), we will limit 
ourselves to present it in subsection |2.2.4[ 


2.2.2 Work of external actions for the considered constrained micromorphic continuum 


According to the procedure shown in the previous section, in order to formulate the equilibrium problem for a 
given continuum, the work of external actions must be given on the portion of the boundary Hj’ = TiTi U ^12 
where tractions are assigned. According to the irreducible expression (19) which we consider for the work 
of internal actions, the external work that balances the boundary terms of the internal one, must contain 
“/orces” / that expend work on the virtual displacement Su and double-forces'" t that expend work on 
the virtual angle variation S(p. In particular, the work of external forces for the considered micromorphic 
continuum is assumed to take the form 


n~)ext _ 

f II — 


{f,5u) ds 


T Sip ds. 


( 20 ) 


Tl 




For the particular case of the bias extension test considered here, we have that traction boundary conditions 
are assigned on the surfaces E3 and E4. More precisely, we have that Eti = = E3IJE4: forces and 

double forces are simultaneously assigned on E 3 and E 4 . More particularly, since in the Bias Extension Test 
E 3 and E 4 are free boundaries the assigned value of forces and double forces is the null value: 


/ = 0 , 


r = 0 . 


2.2.3 Space of configurations for the bias extension test 

Once that the works of external and internal actions have been given for the considered particular case, 
the kinematical boundary conditions must be assigned and the associated space of configurations and of 
admissible variations must be identified. 

For the experimental set-up of the Bias Extension test shown in Fig. we want to assign 

• M = 0 and = 0 on El and 

• u = uq = const and (/? = 0 on E 2 . If we want to frame this situation in the more general picture given 
in section we have to set Ski = ^K 2 = IJ E 2 and the space of configurations Q x D for the 
considered micromorphic continuum subjected to a Bias Extension Test takes the particular form 


Q 

= {u 1 

u = 0 on 

Si 

and 

U = Uq = const. 

D 

= {t 

1 p = 0 on 

Si 

and 

p = 0, on E 2 } 


( 21 ) 


The spaces of admissible variations that must be associated to such space of configurations are given by 

Tu = {Su \ u-\- Su & Q} , Tip = {Sp I p Sp G D} . (22) 

Since both the displacement u and the angle variation p are assigned on Ei IJ E 2 , then it must be <5 m = 0 and 
Sp = 0, so that the spaces of admissible variations can equivalently be written as T„ = {du | = 0 on Ei [J E 2 } 

and Tip = {Sp \ Sp = G on Ei IJ E 2 }. 
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Suitably generalizing what done before, the equilibrium problem for the considered constrained micro- 
morphic continuum subjected to a Bias Extension Test can be formulated as 

,ip*,Su,Sip) + P^^*^{u*,(p*,Su,6ip) = 0 for any {6u,Sip) G 


Find G Q X D such that 


T„ . X T 


ip-, 


where now 7 ?™*, -pext^ rp^ gj,g given by (19l, (20l, (21) and (22) respectively. 

The mathematical question of well-posedness of the geometrically nonlinear micromorphic approach has 
been discussed in ISOl ED EH- the extension of these results to the anisotropic setting is straightforward. 
Some attendant results in the large strain and small strain setting, including an efficient numerical treatment 
and further modeling and well-posedness results can be found in |;19I liHl IHITI IMl EH ES] ■ 


2.2.4 Evaluation of reaction forces and double forces 

As done in subsection 03 we can ask ourselves how, in the framework of the considered constrained 
micromorphic model, we can calculate the reaction forces and double-forces which would be needed to produce 
the kinematical conditions assigned on Ei and E 2 as well as the solution {u*,(p*) obtained solving the 
considered constrained micromorphic equilibrium problem. To do so, we proceed as in subsection 03 
and we introduce a work of external forces also on the portion of the boundary Ei |J E 2 where kinematical 
boundary conditions have been assigned, so that the work of external forces takes the modified form 


Vff = f {f,Su)ds+ f rSipdsF f {fR,Su)ds+ f TaSipds, (23) 

where we remind that = E 3 |J E 4 and that / = 0 and t = 0 for the considered Bias Extension 

Test. Moreover, = E^^ = Ei IJ E 2 . In order to determine the reaction forces and double forces which 
balance the imposed displacements and angle variations on Ei and E 2 , we could pass through the use of the 
irreducible form of the work of internal actions, as we have explicitly shown for the first gradient case in 
subsection 11 . 2.11 

Nevertheless, since for the performed numerical simulations we only used the evaluation of volume inte¬ 
grals, we present here only this last method. To this aim, we write the Principle of Virtual Work for the 
considered particular example, evaluated in the solution {u*, ip*, A*) of the equilibrium problem, as 






Sie - {Kshil - A*) Sis - (a, V(^*) (a, V<5(p) - A* Sip 


- [ - il)SAdv + f {fR,Su)ds+ f TR5ip = Q 

J B 


dv 

(24) 


where the expression for the work of internal actions given in Eq. (18) has been used setting := 14 (u*), 
*6 •= ieiu*), fg := is{u*) and where we impose that this last form of the Principle of Virtual Work must be 
satisfied for any Su, Sip and SA . 

In order to calculate the reaction force, we choose particular test functions Su, Sip and SA such that 


1 . Su is vanishing on Ei non-vanishing and constant on E 2 

2 . Su is an arbitrarily assigned, non-vanishing continuous function outside 'Ski 

3. Sip = 0 everywhere 

4. SA = 0 everywhere. 


With this choice, we can evaluate the reaction force as 


{R, Su) 




Si4 + 



Ste + {Kshil 


A*)Sii 


dv 


(25) 


where we set R := fnds, Si^ := Si 4 (Su), Sig := Sie{Su) and Sis ■= Si 4 {Su). In order to evaluate the 
three components of the reaction force, we choose suitable test functions Su which are aligned with one of the 


17 





















directions of the used reference system. In particular, let 2 3 }) be an orthonormal basis with respect 

to which we want to evaluate the three components Ri = {R, Ci) of the reaction force R. We choose a test 
function Su which possesses all the characteristics previously listed except that the point 1. is replaced by 


6u is equal to on Si . 


In this case the component Ri can be easily calculated according to Eq. (25). 

Analogously, in order to calculate the reaction double-force, we can choose particular test functions Su, 
5(p and 5K such that 


1 . = 0 everywhere 

2 . SCp is vanishing on Ei, non-vanishing and constant on E 2 

3. is an arbitrarily assigned, non-vanishing continuous function outside E/^j 

4. i5A = 0 everywhere. 


With this choice, we can evaluate the reaction double-force as 


TSip = f [—{a,\7p*) {a,\7S(p) — A* Sip] dv, 

Jb 


(26) 


where we set T := Ln Tnds. 

J 2^K2 


2.3 Numerical simulations for the constrained continuum micromorphic model 

In this subsection we show the numerical simulations that we have performed in order to find numerical 
solutions of the equilibrium problem for an unbalanced fabric formulated as the equilibrium of a constrained 
micromorphic continuum subjected to a Bias Extension Test as presented in subsection |2.2[ Such equilibrium 
problem has been implemented in COMSOL® using the weak form package in which the expression (18 1 for 


the work of internal actions can be explicitly given. The work of external actions (20 1 is given in COMSOL® 


just leaving E 3 and E 4 free which means that forces and double forces are assigned to be vanishing on the 
considered subsets of the boundary. Finally, the space of configurations (21 1 is given assigning the kinematical 
boundary conditions on Ei and E 2 as explicitly established in subsection ( 2 . 11 . 


The numerical values used for the constitutive parameters are shown in table ([^. 


Kel 

Ksh 

ai 

CI 2 

0.7 MPa 

21 kPa 

2 kPa-m^ 

0.02 kPa-m^ 


Table 2: Parameters of the constrained micromorphic continuum model. 

Such values have been chosen following a precise heuristic procedure that we present here 

• First of all we remarked the bending stiffness a 2 of the thin yarns is very small (eventually almost 
vanishing), so that we choose a tentative (very small) value for such parameter 

• subsequently, we choose the bending stiffness ai of the thick yarns in order to fit at best the experimental 
S-shape of the specimen (in Fig. |^it is shown as the unbalance on the bending stiffnesses of warp and 
weft is responsible for the S-shape of the specimen) 

• then, we tuned the value of the “sliding” parameter ATei which was seen to have a direct influence on the 
in-plane thickness of the specimen. Turning-on this parameter the height in the middle of the specimen 
itself starts increasing and becomes closer to the real experimental shape (see e.g. Fig. 

• Finally, the in-plane shear parameter Kgh is tuned in order to fit at best the experimental load- 
displacement curve (see Fig. [^. 
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Figure 5: Vertical displacement of the mean horizontal axis for an imposed displacement of 56 mm and 
different values of ai. We remark that no distortion of the specimen is present for a balanced fabric (ai = 02 ). 



Figure 6 : Deformed shape for a displacement of 37 mm and different values of K^i 


We need to spend some extra words on the way in which the reaction force and double force have been 
numerically calculated. We followed the procedure presented in subsection |2.2.4| and, once calculated the 
solution {u*, ip*, A*) for each imposed displacement uq, we calculated the horizontal component of the reaction 
force R with respect to an orthonormal reference frame { 61 , 62 } chosen as in Fig. [^by using equation (25), 
where we chose particular test functions 6u, 5(p and 5A such that 


1. (5u-ei = f, M{X,Y)&B 

2. 6u-e2 = 0, V(V, Y)gB 

5. 6p = 0, y{X,Y)GB 

Y 6A = 0, y{X,Y)€B. 
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Figure 7: Definition of the global ({ 61 , 62 }) and material ({ 7711 , 7712 }) reference frames for the BET. 


The obtained force-displacement curve is shown in Fig. where the comparison with the experimental 
curve and with the reaction force automatically calculated by COMSOL® are also depicted. 



d (mm) 


Figure 8: Load/displacement curve for the constrained micromorphic model 

In figures and [T0| it is shown the deformed shape for two imposed displacements uq, obtained via the 
numerical simulations performed to solve the considered equilibrium of the constrained micromorphic model. 
We remark that the following observations can be easily inferred from the performed numerical simulations: 

• The macroscopic shape of the specimen is completely recovered. 

• The microscopic pattern for the warp and weft is caught to a big extent. Indeed, the real behavior of 
the thick yarns is almost perfectly recovered, while, the thin yarns which are experimentally seen to 
undergo non-negligible slipping, are seen to be fictively elongated. We can affirm that where the thin 
yarns are seen to be elongated with respect to their initial length a slipping is taking place. 

• The shear parameter is chosen to correctly fit the experimental load-displacement curve. We highlight 
the fact that the “force” of the introduced constrained micromorphic model is, to our understanding, 
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the correct one to be compared to the experimental one. In fact, in our micromorphic model, we 
defined “force” the quantity which is dual to the virtual displacement Su. From a direct observation 
of the irreducible form of the work of internal forces (191, we can infer that the internal actions which 
balance a boundary externally applied force are directly determined by equivalent elongations, shear 
angle variations and microstructure-related deformation mechanisms (through the Lagrange multiplier 
A). Analogously, we can infer that the boundary internal actions (double-forces) which are dual to 
the angle variation 6ip are directly determined by the localized bending strains which occur inside the 
considered specimen. We state that the double-forces which are generated by the imposition of the 
angle between warp and weft at the two extremities of the specimen do not contribute to the “force” 
which is measured by the testing machine. Indeed, such angle is maintained fixed by the clamp device 
used during the BET and a supplementary measurement tool would be needed in order to measure such 
sort of couple which is generated by the clamp in order to keep the angle constant at the boundary 
during the test. In summary, we are reasonably assuming that the machine only senses those macro 
and micro deformations modes which induce macroscopic displacements (elongations, angle variations 
and part of the local bending), while the remaining part of the local bending energy which is localized 
at the micro-level is not sensed by the used machine. 



Figure 9: Deformed shape in the simulation (red with black fibers) and experimental (Blue) for a displacement 
of 37 mm 



Figure 10: Deformed shape in the simulation (red with black fibers) and experimental (blue) for a displace¬ 
ment of 56 mm 


As already pointed out, the used micromorphic model possesses a feature that the first gradient models 
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do not posses. In particular, the reaction at the clamps is not limited to a force, but a double-force (a 
couple) arises which is needed to keep the angle between warp and weft constant during the test. Following 
the methods described in subsection 2.2.4 once calculated the solution {u*,(p*,A*) for a given imposed 


displacement uq, it is possible to evaluate the reaction double-force T with respect to an orthonormal reference 
frame { 61 , 62 } chosen as in Fig. [^by using equation (26l, where we chose particular test functions Su, 5(p 
and (5A such that 


1 . ( 5 m -61=0, y{X,Y)eB 

2. 5m • 62 = 0, V(X, Y)gB 

3. y{X,Y)GB 

4. 5A = 0, y{X,Y)€B. 


The results obtained are shown in figure [m where the computed reaction double-force is also compared 
with the one automatically evaluated by COMSOL® . 



d (mm) 


Figure 11: Double-force/displacement curve as reaction T dual of 6(p. 


In a material in which the shear is the main deformation mode one of the most important features to 
check is the angle between the fibers. 
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Figure 12: Angles between the fibers in the constrained micromorphic simulation for a displacement of of 
56 mm 

It can be noticed from Fig. [^that no angle variation occurs in the red regions, while an important and 
almost constant angle variation occurs in the middle of the specimen. The transition from one value of the 
angle to the other one is made thanks to the creation of transition zones in which a smooth variation occurs. 
By direct comparison with Fig. it is possible to remark that the three regions at constant shear angle 
individuated in Fig. are determined by the fact that the thick yarns remain substantially undeformed, 
while the thin yarns are kinematically blocked in the standard triangular zones of the BET, while bend in 
the central part of the specimen. The change of direction of the thin yarns is the only one that determines 
the angle variation between warp a weft. 

We can finally remark that the bending of the thin yarns takes place in well defined transition zones where 
smooth variation of the shear angle occur. 


3 Discrete numerical simulations 


In this section we propose to set up a suitable discrete model in which the motions of the single yarns are 
singularly taken into account. To do so, we decide to use Euler-Bernoulli beams with different bending 
rigidities which are interconnected with rotational and translational elastic springs in order to mimic at best 
the real connections between the yarns. The results obtained with such discrete model allow to 

• explicitly account for the slipping of the yarns 

• better understand the potentialities and limitations of the constrained micromorphic continuum model 
introduced in the previous section. 


The main limitation of the discrete model that we are going to present can be found in the fact that the 
interactions between adjacent fibers are all considered to be elastic. If this can be considered to be reasonable 
to a certain extent, there are for sure some irreversible mechanisms such as friction that cannot be precisely 
described here and that should follow the methods presented in |3n| in order to be fully taken into account. 
Indeed, when unloading the experimental specimen once that the BET has been performed, is not sufficient 
to let the specimen return in its initial configuration. This means that a certain part of the deformation is 
not elastic, but is due to irreversible mechanisms such as friction. Nevertheless a big amount of the imposed 
deformation is recovered and we can hence suppose that the discrete model which we introduce here can be 
thought to be a reasonable compromise between the complexity of the real microstructural motions and the 
simplicity of the model that one wants to introduce. 

The warp and weft yarns are modeled as long Euler-Bernouilli beams disposed at 45 degrees with respect 
to the edges of the specimen (see Fig. 131. The connections between the two families of fibers is guaranteed 
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at the contact points by rotational and translational springs as the ones sketched in Fig. Both beams 
and springs, are considered to be linear. 

More particularly, the two families of yarns are modeled as beams with axial stiffness respectively Ki = 
EAi and K 2 = EA 2 , and bending stiffness respectively = Eli and K 4 = EI 2 where E is the Young 
modulus of the material constituting the yarns, Ai and A 2 are the cross sections and /i, I 2 the equivalent 
moments of inertia. In order to model the interactions between the yarns the set of points in which the two 
families of fibers intersect have been defined in COMSOL®. In this set of points, the interactions between 
the two fibers are supposed to be 

• the resistance to the variation of angle between the fibers, namely the shear stiffness which are accounted 
for by the introduction of a set of rotational springs of stiffness 

• the resistance to the slippages, described by a second set of translational springs of stiffness KgUp to 
describe the response to slippages. 

• the mutual interactions between parallel fibers which are guaranteed by the weaving. Such interactions 
are due to the presence of the orthogonal yarns which connect the two adjacent yarns considered here. 
Therefore a set of springs (modeled as trusses with axial stiffness Kinter) has been inserted between 
every two couples of close interaction points belonging to two different close yarns of the same family 
(see Fig[g. 

A possible downside of such discrete model is its difficult application for bigger specimens due to the high 
number of degrees of freedom needed for a proper description of the yarns response. In this optic, continuum 
models are preferable to discrete ones in view of the design of engineering structures. 



0 0.01 0.02 0.03 0.04 0.05 0.06 0.07 0.08 0.09 0.1 0.11 0.12 0.13 0.14 0.15 0.16 0.17 0.18 


Figure 13: Geometry of the discrete model: undeformed configuration. 

In order to compare the discrete model with the results of the continuous ones it was implemented in 
COMSOL®. In order to have a proper comparison. 

The elastic parameters used are shown in table They have been chosen in order to be reasonably 
compatible with yarns of small cross section area and Young moduli of carbon. The two elongation stiffnesses 
Ki and K 2 have been chosen of the same of magnitude in order to underline the fact that it is indeed the 
difference in the bending stiffness of the two families of yarns which drives the asymmetry of the macroscopic 
deformation. It must be pointed out that as far as this value is high enough it does not have a big influence 
in the results in terms of both displacement and reactions. The parameters relative to the bending stiffness, 
the slipping and the interaction between two fibers of the same set were chosen via a fit of the experimental 
shape of the specimen. In particular the following characteristics were used to fit the different parameters: 
the width of the specimen in the central part, the macroscopic S-deformation, the slipping of the fibers 
and the distance between the fibers of the same set. The shear stiffness was chosen in order to fit the 
experimental force with the reaction evaluated with the simulations. 
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Figure 14: Schematics of the elastic interconnections between warp and weft yarns, (a) rotational spring, (b) 
translational spring, (c) interaction between thin yarns and (d) interaction between thick yarns. 


Ki 

K2 

Ks 

Ki 


K slip 

Kinter 

50000 N 

50000 N 

0.4 N-m^ 

10“^ N-m^ 

2.510-‘‘N-m 

11 N/m 

11 N 


Table 3: Parameters of the discrete model 

As done in the continuum model, in the figures [T5| and [T6| it is possible to see that the discrete model well 
describes the S-response even at different values of displacement. 



Figure 15: Deformed shape in the simulation (green) and experimental (blue) for a displacement of 37 mm 

The fact that the introduced model is able to well describe the slipping between the yarns can be more 
precisely seen with reference to figures and 


25 








































Figure 16: Deformed shape in the simulation (green) and experimental (blue) for a displacement of 56 mm 



Figure 17: Deformed shape for a displacement of 56 mm and some initially superimposed points belonging 
to thick yarns (black) and to thin yarns (red). The highlighted blue points correspond to the experimental 
ones that can be observed in Fig. [T^ 



Figure 18: Deformed experimental shape: individuation of the slipping points. 


In fact, the slippage of thin fibers can be recognized since 


• the points of the thin fibers (red) which are initially located at the free boundaries are finally located 
inside the specimen. 


• the blue points highlighted in Fig. 17 and belonging to the thin yarns are easily recognizable in the 
experimental deformed specimen shown in Fig. and it can be recognized that the simulation well 
describes the experimental behavior, at least qualitatively. 

The response of this model can be compared to the one obtained with the constrained micromorphic 
continuum model. Indeed the different bending stiffnesses of the two families of yarns lead to a set of straight 
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thick fibers and a set of strongly bent thin fibers like in the continuum case. In this discrete simulation, 
nevertheless, the slippages of the yarns are more precisely described. 

We show in Fig. that, with the constitutive choice of the parameters shown in table also the 
load-displacement curve of the discrete model results to be consistent with the experimentally obtained one, 
as well as with that obtained by means of the constrained micromorphic continuum model. We can notice 
that a slightly better fitting is recovered for the discrete model at moderate strains. A better fitting of the 
continuum model xould be obtained for such moderate strains if more complex hyperelastic laws would be 
introduced, but this falls outside the scope of the present paper. 



0 5 10 15 20 25 30 35 40 45 50 55 60 

d (mm) 


Figure 19: Load-displacement curve for the experimental, discrete and second gradient models 

Finally, Fig. shows the angle between warp and weft yarns as obtained with the discrete model., even 
in this model the angle between the set of fibers is almost a constant along the strong fibers direction while 
drastically changes along the think yarns. The results of both the discrete and constrained micromorphic 
continuum models, even with their very different natures, present the same qualitative description of the 
experimental behavior reconfirming the good analogy between the bending stiffness of the yarns and the 
second gradient effects in the continuum. This is another hint toward the importance of the insertion of a 
second gradient energy in a continuous models in order to fully describe the phenomenological mechanical 
response of the woven fabrics. 



A 92.2 
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Figure 20: Angles between the fibers in the discrete simulations for a displacement of of 56 mm 
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Further investigations should be more systematically directed towards the quantitative analysis needed 
to precisely identify the macroscopic constitutive parameters in terms of the microscopic properties of the 
yarns. Suitable multi-scale methods as the one introduced in [5S] may be generalized to be applied to the 
present case. Moreover, the description of the considered system at the microscopic scale may take advantage 
of some of the results proposed in HJ UH |7S1 |70J |771 [751 EHl US] ■ 

4 Conclusions 

In this paper a continuum constrained micromorphic model and a discrete model are introduced to reproduce 
the Bias Extension Test on unbalanced fabrics. We show that both models are able to account for the main 
macroscopic and microscopic deformation mechanisms that take place during the BET up to moderate strains, 
namely 

• the angle variation between warp and weft tows, 

• the different bending stiffnesses of the two families of yarns (which is at the origin of the asymmetric 
macroscopic S-shape of the specimen) 

• the relative slipping of the yarns. 

The results obtained with the two models are satisfactory when considering small and moderate defor¬ 
mations so that it is conceivable to use the proposed models on more extended experimental campaigns. 
This would allow a more precise identification of the introduced constitutive parameters, above all for what 
concerns the different bending stiffnesses which are the main characteristics of fibrous composite interlocks. 

Further studies should be focused on the improvement of the proposed models in the optic of more precisely 
describe irreversible phenomena as friction which can have a non-negligible role during the deformation of 
woven reinforcements. 

Moreover, more extensive experimental campaigns should be carried out in order to determine which is 
the strain threshold until which the integrity of the material is preserved and both the continuum and discrete 
models can be considered to be predictive. In fact, after a given macroscopic deformation, some yarns start 
to be pulled out from the specimen, so that further modeling efforts intrinsically loose their interest. 

Finally, the theoretical and numerical tools presented in this paper can be extended in order to treat 
equilibrium problems concerning 3D composite interlocks. For example, the 3 point bending test on unbal¬ 
anced fabrics could be a useful way of testing and indirectly fitting new constitutive parameters related to 
the out-of plane different bending stiffnesses of the two families of yarns. 
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Appendix A: Alternative numerical implementation of the constrained 
micromorphic model: penalty method 


In this subsection, we briefly mention a method that can be used in order to numerically implement a con¬ 
strained micromorphic model as an alternative to the method of Lagrange multipliers described in subsection 
|2.2.1[ It is known as “penalty method” and consists in implementing a strain energy density which takes the 
form 

W ( 14 , if), ig, (p, ^p') — IV 7 (m 7^65^8) kk/Z ^^coupling (^87 7 

where Wj and Wjj are given in Eqs. and respectively, while the coupling energy takes the form 

K 

^^couplingiis: ~ ^s) 7 


where it! is a constant that may ideally tend to infinity. Indeed, in order to guarantee the boundedness of 
the strain energy density, it follows that p must necessarily tend to ig. We numerically implemented such 
penalty method in order to test the correct convergence of our equilibrium problem formulated with the 
Lagrange multipliers. Being K constant, the considered virtual variations are only 5u and 5p and, moreover, 
the constant K must be chosen sufflciently large in order to guarantee numerical convergence of the solution. 
This last feature can be easily tested by controlling that the solution does not change when increasing the 
value of K (see Fig. 21 1 . The constitutive parameters remain the same as the ones used in the numerical 
simulation with the Lagrange multiplier (see table|^. It is possible to notice that, suitably increasing the value 
of K the model converges and the limit corresponds to the solution obtained with constrained micromorphic 
simulation with Lagrange multipliers (see also Fig. |^. 



Figure 21: Vertical displacement of the mean axis for a displacement of 56 mm and different values of K 

This penalty method can be seen as a useful tool for the easy implementation of constrained micromorphic 
models due to their high numerical stability. 
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